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Abstract. Quantum control is traditionally expressed through bilinear models 
and their associated Lie algebra controllability criteria. But, the first order 
approximation are not always sufficient and higher order developpements are 
used in recent works. Motivated by these applications, we give in this paper a 
criterion that applies to situations where the evolution operator is expressed as 
sum of possibly non-linear real functionals of the SEime control that multiplies 
some time independent (coupling) operators. 
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1. Background on quantum control 

Controlling the evolution of molecular systems at quantum level has been envi- 
sioned from the very beginings of the laser technology. However, approaches based 
on designing laser pulses based on intuition alone did not succed in general sit- 
uations due to the very complex interactions that are at work between the laser 
and the molecules to be controlled, which results e.g., in the redistribution of the 
incoming laser energy to the whole molecule. Even if this circumstance initially 
slowed down investigations in this area, the realization that this inconvenient can 
be recast and attacked with the tools of (optimal) control theory [TS] greatly con- 
tributed to the first positive experimental results [U [5T1 [331 El 13 EZl HO] ■ 

The regime that is relevant for this work is related to time scales of the order 
of the femtosecond (10^^^) up to picoseconds (10^^^) and the space scales from 
the size of one or two atoms to large polyatomic molecules. 

Historically, the first applications that were envisionned were the manipula- 
tion of chemical bonds (e.g., selective dissociation) or isotopic separation. Although 
initially, only few atoms molecules were investigated (di-atomics) the experiments 
soon were designed to treat more complex situations [2 as selective bond disso- 
ciation in an organi-metalic complex CpFe{C0)2Cl {Cp is the cyclopentadienyl 
ion) by maximizing or minimizing the quotient of CpFeCOC'l^ ions obtained with 
respect to FeCl^ ions. 

Continuing this breakthrough, other poly-atomic molecules were considered 
in strong fields. For instance, in [2T] the molecules are the acetone {CH-i)2C0, the 
trifluoroacetone CH^COCF^ and the acetophenone CeH^COCH^. Using tailored 
laser pulses it was shown possible to obtain CH3CO from {CH^)2CO, CF3 (or 
CH3) from CH3COCF3 but also CeH^CH^ (toluene) from CeH^COCH^. 

But the applications of laser control do not stop here. High Harmonic Gen- 
eration) [7] is a technique that allows to obtain output lasers whose frequency is 
large interger multiples of the input pulses. 

A different class of applications works in a different regime of shorter time 
scales and large intensity. This regime is additionally not compatible with the 
standard Born-Oppenheimer approximation and requires to consider both nucleari 
and electrons as quantum particles with entangled wevefunction [IJ . 

In a different framework, the manipulation of quantum states of atoms and 
molecules allows to envision the construction of quantum computers [T31 [57] 

Finally, biologically related applications are also the object of ongoing re- 
search. 



2. Background on controllability criteria 

We start in this section to investigate the theoretical controllability results that 
are nowadays available for quantum systems. The evolution of the system will be 
described by the driving Schrodinger equation (we work here in atomic units i.e. 
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h=i) 

i^^^{t,x)^H{t)^{t,x) (2.1) 

where H{t) is the Hamiltonian of the system and x G the set of internal degrees 
of freedom. We introduce the Hilbert space structure given by the scalar product 



if, 9)= / f{x)g{x)dx (2.2) 



where a + ib = a — ib the conjugate of a complex number. 

We only consider in this paper situations when the Hamiltonian is auto- 
adjoint H{ty = H{t); we denoted by T^' the adjoint of a operator T. The auto- 
adjointeness of H implies that the L'^{'EO) norm of the evolving state is conserved. 
Indeed 

^||*(X,<)|U2(R.) = j^{^{x,t),^{x,t)) 

= {j^^ix, t), ^{x, t)) + {^{x, t), j^^ix, t)) 

= {Wiif{x,t),^{x,t)) + {^{x,t), ^^{x,t)) = 0. (2.3) 
t I 

Thus 

||*(a;,i)||i,2(R7) = ||*o|U2(RT), > 0, (2.4) 
so the wave function ^'(t), evolves on the (complex) unit sphere 

S={i,&L\W) : ||^|L.(R.) = l}. 

When the system evolves freely under its own internal dynamics i.e. when 
isolated molecules are considered, the free evolution Hamiltonian Hq is introduced. 
This Hamiltonian is the sum of the kinetic part T and the potential operator V{x) 
: Ho = T + V{x). A prototypical example of T is the Laplace operator while for 
V{x) one can encounter Coulomb potential or Lennard- Jones type dependence. 
We obtain the following evolution in the absence of external interaction: 

i^^^{t,x) = Ho^{t,x) (2.5) 
^{to,x) = ^o{x). 

But, when the free evolution of the system docs not generate a satisfactory dynam- 
ical output, an external interaction is introduced to control it. An example of exter- 
nal control of paramount importance is a laser source of intensity e{t) G R, t >0. 

The purpose of control may be formulated as to drive the system from its 
initial state \E'o to take a convenient dynamical path to a final state compatible 
with predefined requirements. The control is here the laser intensity e{t). We will 
come back later with details on the laser field e{t). 
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This laser will modify the Hamiltonian H{t) of the system. A first order ap- 
proximation can be considered by introducing a time-independent dipole moment 
operator fJ.{x) resulting in the dynamics: 

i^^^it,x)^{Ho + e{t)^l)^it,x) (2.6) 
*(io,2;) = "foix). 

This is the so-called bi-linear framework (the control enters linearly multiplying 
the state), that is the object of much theoretical and numerical work in quan- 
tum control. We also review below some of the results that are available in this 
formulation. However, recently, higher order field dependence has been consid- 
ered in different circumstances see e.g., [THITn] for details. In these situations the 
Hamiltonian H{t) is developped further as : 

H{t) = Ho + e{t)ni + e2(i)^2 + ... + e^{t)^L- (2.7) 

The question that will be of interest to us in this work is the study of all 
possible final states for the quantum system. This question is important in order 
to understand the capabilities that a laboratory experiment will be able to pro- 
vide and also, in a more general setting, to accompany the introduction of new 
experimental protocols. 

More specifically, we will show how the criteria available for bilinear control 
can be extended to treat the Hamiltonian (|2.7p where a single control amplitude 
e{t) appears before different coupling operators /Lti, 

Many of the questions regarding the properties of the quantum control pro- 
cedures, such as controllability, optimal control definition, etc, ... need, in order to 
be defined, to specify the admissible control class, i.e., the set U where the control 
e(t) is alowed to vary. Among the properties that can define this admissible set, 
some are related to the regularity of the time-dependence (L^, 7?^, ... etc) or of the 
Fourier expression (sum of sinusoidal functions multiplied by an overall enveloppe, 
etc,...) or to additional structure: e.g. piecewise continuous, piecewise constant, 
locally bounded ... 

The choice of one or several conditions in the list above is motivated in 
practice by capability to reproduce that particular form or to inherent experi- 
mental restrictions (finite total laser energy /fluence, etc). As the laser technol- 
ogy is constantly evolving, the first class of constraints becomes less critical and 
thus it is realistic to consider very weak constraints on the control set, e.g. U = 

i2(R)ni-(R). 

However, to treat even more general situations, we will consider in this work 
controls e{t) that are piecewise constant, taking any value in a set V ^ which will 
remain fully general. 

2.1. Infinite dimensional bilinear control 

When compared to the finite dimensional control equations (see Section [^^ . con- 
trollability of the infinite dimensional version of the bilinear Time Dependent 
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Schrodinger Equation is much less understood at this time. In fact, most of the 
progress obtained so far takes the form of negative results, in contradiction with 
the positive results available in finite dimensional settings. However we see the ab- 
sence of positive controllability results is rather a failure of today's control theory 
tools to provide insight into controllability rather than an actual restriction. We 
do believe that new tools and concepts will make positive results possible. 
Let us write the solution of (|2.6p in the following form: 



This formulation (see [TT] for details) is granted by the properties of the operator 
/i : Hq{IV) — > H~^{IV) which is continuous when /i is bounded; we also recall 
that the control e can be considered bounded in both L°° and L^. 

The application e(t) vl/(x,t) possesses an important compacity property 
which is the key of the controllability results (we refer the reader interested in 
details to [1[25]) : 

Lemma 2.1. Suppose that ii : X ^ X is a bounded operator and that Hq generates 
a C° semigroup of bounded linear operators on some Banach space X (e.g. X = 
H^ilC)). Denote for T > Q and e G L^([0,r]) by '^,{x,t) the solution of ([^ 
with control e. Then e "^f^ is a compact mapping in the sense that for any e„ 
that converges weakly to e in L^{[0,T]) '^^^^ converges strongly in C{[0,T]; X) to 



This compactness property allows to give negative results for general bilinear 
controllability settings as in [5] where they were applied to the wave and rod equa- 
tions. Specific statements for quantum control have been latter derived (Thm. 1 
from [5D] ; see also [31 [52]) and can be stated as: 

Theorem 2.2. Let S be the complex unit sphere of L^{K^). Let ^ be a bounded 
operator from the Sobolev space X ( e.g., X — {IV) ) to itself and let Hq generate 
a semigroup of bounded linear operators on X. Denote by 5'e(2;,t) the solution 
of (|2.6p . Then the set of attainable states from ^I'o defined by 



is contained in a countable union of compact subsets of X . In particular its com- 
plement S n X \ AS with respect to S Ci X is everywhere dense on S Ci X . The 
same holds true for the complement with respect to S . 

In a different formulation, the theorem implies that for any ^E'o G X fl 5, 
within any open set around an arbitrary point \1/ S X n 5 there exists a state 
unreachable from '^q with L^ controls. 

Remark 2.3. Note that the result does not give information on the closure of the 
set ^5*. In particular it may well be that while AS still has dense complement its 
closure be the whole space X. This would be the so-called approximate controlla- 
bility i.e. the possibility to reach targets arbitrarily close to any given final state. 




(2.8) 



AS = Ut>o{*.(x, T); e(t) e L\[0, T])} 



(2.9) 
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Despite some attempts in the litterature, at this time there is no answer (positive 
or negative) to this question. Among the ingredients that make this study difhcult 
we can mention the possibihty to use arbitrary large final time T, the necessity to 
treat the continuous spectrum of the operator Hq and the intrinsically unbounded 
domain on which the problem is posed. 

To complicate even more the landscape, situations exists where the results 
obtained in infinite and finite dimensional representation are of different nature. 
We will illustrate with a classical result on the harmonic oscillator. 

Lemma 2.4. The infinite dimensional harmonic oscillator Hq = —-§^ + x , fj. = x 
is not controllable. Moreover the set of all admissible states is a low- dimensional 
manifold of . 

Proof. Let us begin by noting that the operators —iHo and —ifi form a Lie algebra 
of dimension 4. Indeed, let us compute the iterated commutators of Ho = —-§^+x 
and jjL ~ X : 



Thus the dimension of the Lie algebra — 4 and as such the system cannot be 
controllable (all the states are on a low dimensional manifold of L?). We refer 
to f23^ for recent contributions when the algebra of the operators Hq and /i is 



This result is to be contrasted with additional works that show that any 
(spectral) truncation of the harmonic oscillator is controllable (see [26] for details) . 

What can be deduced from the above result is that truncating an infinite 
dimensional system is not always justified and care must be taken to check that 
the control obtained in the resulting finite dimensional approximation remain a 
good control for the initial, infinite dimensional system. Of course, this is not 
needed for situations which are inherently finite dimensional quantum systems 
(e.g., spins). 

2.2. Finite dimensional bilinear control 

Here our focus will be on finite dimensional systems. We introduce an orthonor- 
mal basis D = {'il)i(x)]i = 1, ..,-/V} for a finite dimensional space. An important 
example of such a space is the one spanned by the first N eigenstates of the in- 
ternal Hamiltonian Hq. This example is also motivated in bi-linear settings by 
the "perturbation" argument that considers the control term e(t)/x as a first order 
developpement of H{t) around Hq. Note however that no concept of "smallness" 
is introduced in the definition of admissible controls lA. 




(2.10) 



(2.11) 



(2.12) 



finite dimensional. 



□ 
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Denote by M the linear space that D generates, and let Ho-a.b ~ {Hotpa, V'fc) 
and fie-a.b = {'>jJa, fJ't'tl'b) be the expressions of the operators Hq and fi£ with respect 
to this basis, £ — 1, ...,L. To keep notations simple we will still denote from now 
on by Hq and fi{ the resulting N x N symmetric matrices. 

In the Galerkin approach, expressing the Schrodinger equation in the space 
M is equivalent to supposing '^{x,t) = X^i^^i 4'i{^)ci{t). 

^ dc{t,e,co) ^ ^_ ^ ^^^^^^^ ^ ^ e^{t)fXL] c{t- e; cq) (2.13) 

c{t = O;e;co) = cq. 

In the following, when no ambiguity prevents it, we will also simply denote c — 
c{t; e; cq). The finite dimensional counterpart of the norm conservation property (12. 4p 
reads: 

N 

i.e., the state c evolves on the unit sphere Sn of C'^. The controllability can be 
formulated in this case as: 

Definition 2.5. The system {Hq, /ii, /il) is called (wavefunction) controllable, if 
for any two states Ck G Sn, fc = 1, 2 there exists a final time T < oo and control 
e{t) € I/^([0,T]) such that the solution of eqn. (|2.13|) starting from ci ends in C2 
at final time T: c{T] e; ci) — C2. 



Although specific results for this setting exist [3T1 [32] , a different alternative 
is to see (|2.13p as a system posed on [/(A^ 0- We introduce the evolution equation 
on U{N): 

'^^^dt^ = [Ho + eW/ii + ... + e''{t)^iL] U{t; e) (2.14) 
U{t = 0;e) = Id. 

Since Hq and fii are symmetric matrices, U{t; e) will remain unitary for all t > 0. 
It is classical to remark then that the evolution of c{t; e; cq) can be obtained from 
the evolution of U{t; e) by 

c(t;e;co) = U{t;e)co. 

In particular it follows that if the set of all attainable matrices U{t; e) is at least 
SU{N) then the system is controllable. This is almost a necessary condition for 
controllability, a notable exception being the circumstance when N is even: in this 
case, if the set of all attainable matrices contains Sp{N/2) then controllability still 
holds. We refer to [101 1^1^ for more detailed information. 



U{N) is the set of all Af X A'' complex unitary matrices. 
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Let us just mention that different representations of the system include the 
density matrix formulation with time dependent density matrix operator p{t) sat- 
isfying 

d 

z— e; po) = [i?o + e{t)pi + ... + {t)pL, p{t; e; po)] (2.15) 
p{t = 0; e; po) = po 

Then one can show p{t]€;po) — U [t] e) p^W [t; e) . Controllability in this case is 
the possibility to steer any initial mixed state po to any other state p^ unitarily 
equivalent to 110. 

Note that the density matrix controllability is equivalent to requiring that 
the set of all matrices attainable from identity be at least SU (N) . 

At a general level, the evolution equation (|2.14|) can be re-written as 

^^IM^ ^^A + e{t)B, + ...e\t)BL)x{t- e- xo) (2.16) 

a;(0)=a;o. (2.17) 

where x{t;e;xQ) belongs to a Lie group G ( see [16j |8j [9] for basic facts about 
the lie groups) and A,Bi,...,Bk to its associated Lie algebra L{G). The equation 
above is to be taken in the usual sense (using the exponential map) when e.g., e{t) 
is piecewise continuous/constant and in a weak sense (integral form) for general 
e(t) (see e.g. [3; for additional details). For the quantum control problem A = —iHo 
and Bi = -ipe, G = U{N). 

Remark 2.6. Everything that will be said in this and following sections applies 
with trivial modifications to the situation of several laser fields. For notational 
convenience we will only give here the results for a unique laser field. 

We will denote by LA,Bi,...,Bk C L{G) the Lie algebra spanned by A, Bk, 
k — 1, K and by e the unity of G. 

Let us now consider the set of all reachable states from an initial state y: 

T^uiv) = {2;(i; e; V) solution of (l2J6l) ; e £ W}. (2.18) 
It is immediate to see that 

K{v)^nii{e)y (2.19) 

and thus, describing the set 7^^(e) allows to completely describe all other reachable 
sets. When the final time is not specified, we will denote 

nuiv) = ^t>oK{v)- (2.20) 

The central question is to characterize TZjj{e). When the bi- linear setting is con- 
sidered i.e. L = 1 and we note B — Bi, we have the following result [TOl 



TV X TV matrix p2 is said unitarily equivalent to a TV X TV matrix pi if there exists M G U(N) 
such that p2 = TUpoAft. 
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Theorem 2.7. Consider the system (j2.16p defined on a Lie group G with associated 
Lie algebra L{G) containing A and B . If G is compact and the Lie algebra La.b 
generated by A and B is the complete algebra L{G) : La,b = L{G) then the set 
TZ^{e) of all states states from the identity is the Lie group G. Moreover, there 
exists < T < oo such that TZj^ (e) = G for all T' >T . 

This gives, when applied to quantum control [2A : (L = 1, /i = /ii): 

Theorem 2.8. // the Lie algebra L^iH„.~-ip, generated by —iHo and —ifi has di- 
mension N"^ (as a vector space over the real numbers) then the system Ij2.14p is 
density matrix controllable. Furthermore, if both ^iHo and —ifi are traceless then 
a sufficient condition for the density matrix ( thus wavefunction ) controllability of 
quantum system is that the Lie algebra L^iH„,-i^ has dimension N"^ — 1. 

Although the results above conveniently address the situation of a bi-linear 
setting, we are not aware of any similar results for the general quantum control 
situations (|2.7p . In particular, we know by the result above that, if ui,...,ul are 
independent controls, i.e., 

'^''^^'^^p^ ^ (A + u^{t)B, + ... + UL{t)BL)x{t- e; x^) (2.21) 
x(0) = xo, (2.22) 

an equivalent condition for the controllability of the above system on its compact 
Lie group G is that A,Bi,...,B]^ generate the whole Lie algebra L{G). But, there is 
no obvious way to say what will happen when the controls ui are not independent 
but related by the condition ue = e^(t). This study is the purpose of the next 
section. 



3. Criteria for non linear operators 

In order to extend the controlability results above beyond bi-linear interaction 
Hamiltonians, we will introduce in this section a more general setting: we will 
rewrite the control equation (I2.16P as 

^^^^^pl = {F,{eit))Bi + ... + Fi{eit))BL)x{t; e; Xo) (3.1) 
x{0) = xq. (3.2) 

where Fk : V M, are real functionals. Note in particular that we do not impose 
any assumption on the regularity of the functionals F^. Of course, one can recover 
the equation (|2.16p by setting Fk{x) — and adding F^ — \. 

In order to avoid trivialities, we will suppose in the following that 

the functionals (-Ffc)^^! are linearly independent. (3.3) 

Otherwise one may just consider a subset of functionals that are linearly indepen- 
dent and adjust the matrices B^. accordingly. Of course, since we do not specify 
the set V that lists all the possible control values e the hypothesis above needs to 
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be understood in the foUowing acccption: the functionals Fk are said to be hnearly 
dependent if there exist constants Ai, £ R such that J2f=i "^j-^ji"") ~ ^'^^ 
all V G V. Otherwise the functionals are said to be linearly independent. 

In order to obtain the quantum controllability results, we begin in this section 
with a controllability criterion on compact Lie groups. These results build on 
classical references for bilinear controllability [191 . We give first a weak but intuitive 
form and then we state the fully general one. 

Theorem 3.1. Let (j3.ip be a control system posed on a compact connected Lie 
group G, with linearly independent functionals (-Ffe)^^x- Then if the Lie algebra 
generated by Bi,...,Bl is the full Lie algebra L{G) of the group G, then the system 
is approximatelly controllable, i.e. for any a,b G G, b is an accumulation point of 
the set of all states x{t) attainable from x{0) — a with admissible controls. 

Proof. Let us begin by noting that if Fk are independent then there exist values 
ej € V, j = 1,...,L such that the vectors v{ej) — (i^i(ej), ^^(ej)) are hn- 
early independent. Suppose on the contrary that this is not true. Consider then 
a maximal set of vectors v{Ei), ...,v{Ep) that are linearly independent. The ma- 
trix iFk{Ej))j^!^i.j^i has rank preciselly p and thus one can extract p functionals, 
denoted for notational convenience Fi, Fp such that 'rank{Fk{Ej))^'^^.j^^ = p. 
Take now some functional Fp+i not in this set. It follows that rank{Fk{Ej))''^'t\]^^\ = 
p and as such det{Fk{Ej))^'t\'^^\ = for any Ep+i G V. This determinant can be 
computed as: 

det{Fk{E,))ltiPtl = XiFi{Ep+i) + ... + Xp+iFp+i{Ep+i) = 0. (3.4) 

Note that A^ do not depend on i?p+i and that in particular 

Xp+i^det{Fk{E,))';:t,.^^^,^0. 

Thus eq. p.4p implies that a linear combination with at least one non-null coeffi- 
cient Ap+i exists such that XkFk{E) — for all E €V. This is prevented by 
hypothesis. 

We have thus proved the existence of Cj G j = 1, ...,L with the v{ej) = 
(Fi(ej), i^L(ej)) linearly independent. This means that Mj = Ylik=i^k{ej)Bk 
are also linearly independent and span the same linear space as Bk, k — 1,...,L 
and thus Mj span also the Lie algebra L{G). Moreover, all states {e**^Ja;(0); t £ 
j < L} are atainable from x(0) for the control system (|3.ip . 

It is clear that to prove approximate controllability is sufhcient to set a — e 
the neutral element of the group G. i.e. we have to prove that the closure TLu (e) 
(with respect to the Lie group topology) of the reachable states from identity is 
the whole G. From the hypothesis and sujectivity of the exponential mapping this 
is equivalent to proving that 



{e^'-M €L{G)} cnu{e). 
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We will begin by noting that TZu(e) is a group. Indeed, take two elements 
x{ti;ei;e),x{t2;e2,e) E TZu[e). Then,defining the control ei2 : [0,ti + ^ R by 
£12 (i) — ei(i) for all < t < ti and £12(^1 +t) = e2{t) for all < t < t2 we obtain 
x{ti+t2; e; e) = x{t2', €2; e)x{ti; ei, e) and thus x{t2', £2; e)x{ti; £1, e) € TLyie). Hence 
'R-y{e) is a semi-group which implies that TZuie) is a semi-group too. 

Let us now consider a G TZu{e). Then a" S TZuifi) for any n = 1,2,.... 
Since TZu{e) C G which is a compact group, TZu{e) is compact at its turn. Then 
there exists a sequence, that we can take such that Uk with Uk — Uk-i > 2, with 
a"fc ^ ft £ G. But then ^J7(e) 3 a"'=""'=-i"i ^ bb-^Q-^ and thus a^^ £ ^57(e). 

It is immediate to see that, since the solution for the control e{t) = Cj is 
x{t\ 0, e) = e**^^ we have the inclusion {e**^^ ; < > 0, j < L} C nijfi). Since 
is a group, we will also have {e**^^ ;t € R;j < L} C TZu{e). Consider now two 
matrices Xi,X2 G i(G) such that 

{e*^';i > 0} C TZu{e), i = 1,2. 
We invoke now the formula 

^t[Xi,X2] ^ /'g-^X2/V^g-^Xl/v^gtX2/v^g^^l/^A^^"■ ^3 5^ 

to conclude that 

Similarly, we use the formula eJ^^^+*'^^^ = lim„^oo (e'l-^i/^e*^-^^/")" to con- 
clude that 

{e'''''+''''';h,t2 e-R} cniA^. 

We have thus proved that the set {M e L{G);e*^^ £ Tiu{e); £ R} 
contains Mj, j = 1,...,L, is closed to commutation and is a real vector space. 
Thus it contains L{G) hence the conclusion of the theorem. □ 



The Theorem above has the advantage to be both intuitive and self-contained. 
However it only gives approximate controllability results, which are not the strongest 
forms available. But, in order to obtain exact controllability more involved tech- 
niques are needed. In the litterature, similar situations are treated by making use 
of the Chow theorem [12] and of the bi-linear control techniques [HI [28]. The 
criterion can be stated as follows: 

Theorem 3.2. Let ()3.H) be a control system posed on a compact connected Lie group 
G with linearly independent functionals Fk : V — > R, k = l,...,L and piecewise 
constant controls e taking any value in some set V . Then a necessary and sufficient 
condition for the exact controllability is that the Lie algebra Lbi,...,Bl generated 
by Bi,...,Bl be the full Lie algebra L{G) of the group G. 
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Proof. We recall (see also end of Section I2.2p that the set of attainable states is 
included in the set of attainable states for the system 

'^"'^^i^';''"^ = [uiit)B^ + ... + UL{t)BL]xit; e; xq) (3.6) 

x{0)=e, (3.7) 

whose controllability is equivalent to ^^Lbi....,Bl = L{Gy\ Thus Lbi,...,Bl ~ L{G) 
is a necessary condition for controlability. To prove that is also sufficient, consider 
as in the proof of the Theorem lS.ll the matrices Mj = Fk{e.j)Bk, j = 1, L 

that generate the same Lie algebra Lbi,...,Bl- We recall that all states {e**^j;i G 
]R+, J < L} and all finite products of such states are atainable from the identity e. 
We invoke now a technique present in the proof of Thm. 3.1 of [28]: for any P G IN 
and any multi-index i — (ii,...,v) G {1,...,L}'' denote by A{i,T) the atainable 
states with the sequence of operators i and total time less than T : 

{r r ^ 

e=i 1=1 ) 

We know by the Chow theorem that the union of the sets A(i, T) is the whole Lie 
group G. Also, it is immediate that any A{i,T) is image of a compact set thus 
compact. It follows by the Baire category theorem that A{i,P) has non-empty 
interior at least for a couple {i,P). For such an i — (ii, ...,?to) we introduce the 
mapping P : R™ G defined by t ^ {ti,...,U) ^ F{t) = YW^^e^'^'i . This 
mapping is analytic and its image is has nonempty interior. By the Sard theorem 
its differential dF{t) has full rank (i.e. equals the dimension of the tangent space 
TG of G) at least at some point t and thus in a neighborhood. But since dF{t) 
depends analytically on t the set of points where the rank is full is dense in R™ 
and as such the rank is full for some t with all components strictly positive. Using 
a local inverse mapping theorem it follows that the image F{T) has non-empty 
interior where T is an open subset of R™. But all points in F{T) are realisable 
with admissible controls and thus the set of reachable points TZjj{e) contains an 
open subset D of G. 

By the previous Theorem, 7^^(e) is a subgroup i.e. for any y G Ti.i({e) the 
set Dy is also reachable. Since in addition Tli^{e) is dense in G it follows that 
n^{e) = G. □ 



4. Applications to quantum control 

The purpose of this section is to instantiate the results obtained previously to the 
specific situation of the quantum control. We will give two results, one for the 
density matrix formalism and the second for the wave function. 

4.1. Density matrix 

To consider the specific situation of the density matrix formalism, we use the 
results of the Section [3] for the Lie group U{N). We obtain a first 
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Theorem 4.1. Consider the system 

i^^p{t- e; po) = [^^0 + FMt))^^l + - + FL{e{t))^iL, p{t- e; po)] (4.1) 
p{t = O;e;po) = Po 

and suppose that the family {1, i^i, i^^} is linearly independent. 

Then, when at least one matrix Ho,fJ,i,---,^iL has nonzero trace, the equa- 
tion (|4.ip is density matrix controllable if and only if the Lie algebra LiHoAfi^ ifi^ 

spanned by the matrices iHo,ipi,...,ifXL 'is the Lie algebra u{N) of all skew-hermitian 
matrices or equivalently dimYiLiHa,i^i,...,ii_LL = 

Otherwise, when all matrices Hq,ij,i,...,pl have zero trace, a necessary and 
sufficient condition for controllability is that LiH„,ip,i,...,ip,i^ — su{N) or equivalently 
dimjiLiH„,t^i,...,t^L ^ N'^ -1. 

Proof. The first part of the conclusion follows from Theorem 13. 21 for the Lie group 
G^U{N). 

When all matrices have zero trace one uses the same result for G — SU (N) 
noting that if two matrices pi and p2 are unitarily equivalent p2 = MpiM^ then 
there exists ^ e R with AI^u = Me'f e SU{N) and p2 = MsuPiMj^. □ 

An algorithmic verification of the above theorem can be devised as follows: 

1. Test whether the functions {1, Fi, F^} are linearly independent. If the 
answer is yes go to next step, otherwise keep only a subset Fi-^,...,Fi^ with 
{1, Fi-^, Fi^} linearly independent and modify the Bi ... ,Bl accordingly. 
For notational convenience we suppose all functionals are independent (p — 
L). 

2. Construct the traceless matrices Hq = Hq — "^'^"'' /d, pi = fii — ^^^^p^Id,..., 

PL ^ pi - ^''j^^^ /d. Denote by O = {iHf),i'p{, ...^ipi] . 

3. Write any element of O as a column vector and compute the rank r — 
rank{0) over the real numbers. 

4. Construct all commutators C of matrices in and test whether rank{0 U 
C) = r. If not, set O := O U C and return to previous step. 

5. Test whether r = iV^ — 1. If yes the system is controllable, if not the control- 
lability does not hold. 

Even more precise results can be derived for the situation in Eqn l2.15l 

Theorem 4.2. Consider the developpement of the interaction Hamiltonian H — 
Hq + e{t)pi + ... + e^{t)pL resulting in the following evolution equation 

d 

i—pit; e; po) = [Ho + t{t)pi + ... + e^{t)pL, p{t; e; po)] (4.2) 
p{t = O;e;po) = po 



^Some optimizations are possible at this point as only new commutators are generally needed to 
be computed. We do not enter into details here. 
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Then, when at least one matrix Hq,^i,...,pll has nonzero trace, the equa- 
tion (14. 2p is density matrix controllable if and only if the Lie algebra LiHo,ip,i,...,ip,,^ 
spanned by the matrices iHQ,i^i,...,i^L is the Lie algebra u{N) of all skew-hermitian 
matrices or equivalently dimnLiHo.ini....,inL — 

Otherwise, when all matrices -ffo^/^iwML have zero trace, a necessary and 
sufficient condition for controllability is that Lifj^j i^-^^ i^^ — su{N) or equivalently 
dimnLiHo,itJ.i,...,itJ.L = iV^ - 1. 

4.2. Wave function 

To derive results for the wave function of the same nature as the two criterions 
above one has to analyse the transitive subsets of U{N). We recall that a subset 
A C U{N) is called transitive when for any two vectors a, b on the unit sphere 
of C ^ there exists a matrix X G A with b — Xa. For the situation of quantum 
control, such a study is available in the literature [1]. To be able to state the 
corresponding result for this specific situation here, we introduce the centralizer 
Cqz of an element z ^ G which is defined as the set of all elements that commute 
with z: 

Cgz = {x E G : xz = zx}. 
We also define P^i- diag{l, 0, 0) E U{N). 

Theorem 4.3. Consider the system 

^dc^t^epl ^ ^ FMt)hi + - + FLm)l^L] c{t; e; cp) (4.3) 
c{t = O;e;co) = cq. 

with \\cq\\ — 1. Suppose that the family {1, Fi, F^} is linearly independent and 
denote by ii_f/o,iMi> - >vt algebra spanned by the matrices iHo,ifj,i,...,i^L. 

Then the equation (j4.3p is ( wave function ) controllable if and only if 

dimYiLiHo,iiii,...,iiiL ~ dim{LiHo,ifj.i,...,ifj.L ^ CqP) = 2iV - 2. 
In particular a sufficient condition for controllability is that 

dimY(LiHoAtii,...,ifiL = 
Proof. The proof follows from arguments in [T]. □ 

The following procedure allows to implement the above criteria: 

1. Test whether the functions {1,^1,...,^^} are linearly independent. If the 
answer is yes go to next step, otherwise keep only a subset Fi^,...,Fi^ with 
{1, Fij^ , Fip} linearly independent and modify the Bi ... ,Bl accordingly. 
For notational convenience we suppose all functionals are independent i.e. 

2. Denote O = {i-ffo, ^Mii ■••I ■ 

3. Write any element of O as a column vector and compute the rank r — 
rank{0) over the real numbers. 
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4. Construct all commutators C of matrices in c|f| and test whether rank{0 U 
C) = r. If not, set O := O U C and return to previous step. 

5. Extract from O the matrices that commute with P and compute the rank d 
of this ensemble over JR. Test whether r — d = 2N — 2. If yes the system is 
controllable, if not the controllability does not hold. 

We also obtain 

Theorem 4.4. Consider the developpement of the interaction Hamiltonian H = 
Hq + + •■■ + resulting in the following evolution equation 

^dcit^epl - [ffo + <t)^i^ + ... + e^WML] c{t- e; cq) (4.4) 
c{t = 0;e;co) = Cq. 

with \\cq\\ — 1. Denote by ii_ffo,ipi,...,ipi, the Lie algebra spanned by the matrices 
iHo,ifJ,i,...,ifJ.L- 

Then the equation (j4.4p is ( wave function ) controllable if and only if 

dimuLiHo,iiii,...,iiiL - dim{LiH„Atj.i,...Atj,L f^CaP) ^2N -2. 
In particular a sufficient condition for controllability is that 

dimfiLiHo,ini,...,inL = N"^ ■ 

Remark 4.5. All the above results basically state that controllability with linearly 
independent functionals of a single control e is true whenever the same equation, 
but with completelly independent controls, is controllable. 
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